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Abstract. We present an empirical determination of the mass-loss rate as a function of stellar luminosity and 
I | effective temperature, for oxygen-rich dust-enshrouded Asymptotic Giant Branch stars and red supergiants. Fo 

this aim we obtained optical spectra of a sample of dust-enshrouded red giants in the Farge Magellanic Cloud, 
' which we complemented with spectroscopic and infrared photometric data from the literature. Fwo of these 

, turned out to be hot emission- line stars, of which one is a definite B[e] star. Fhe mass-loss rates were measured 

1 through modelling of the spectral energy distributions. We thus obtain the mass-loss rate formula log M = 

—5.65 + 1.05 log(L/10, 000 Lq) — 6.3 log(T c ft/3500 K), valid for dust-enshrouded red supergiants and oxygen-rich 
AGB stars. Despite the low metallicity of the LMC, both AGB stars and red supergiants are found at late spectral 
in . types. A comparison with galactic AGB stars and red supergiants shows excellent agreement between the mass- 

loss rate as predicted by our formula and that derived from the 60 /im flux density for dust-enshrouded objects, 
but not for optically bright objects. We discuss the possible implications of this for the mass-loss mechanism. 

Key words, stars: AGB and post-AGB - stars: carbon - stars: mass-loss - supergiants - Magellanic Clouds - 
Infrared: stars 
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1. Introduction tive envelope, but episodically the temperature becomes 

sufficiently high to temporarily ignite a helium-burning 

Stellar death results in the injection of chemically en- gheU (thermal pulges) This leadg tQ penetration of the 

< , riched mass, energy and momentum into the interstellar ,. , n r i j j. 

J3 , convective envelope, allowing tor nuclear processed mate- 

CO , medium. It also results in baryonic mass being locked up ^ tQ bg transported to the steUar surface (e g ? carbon 

forever into degenerate objects (white dwarves, neutron , , , , . jit,- \ 

, , , . ' ana s-process elements sucn as lithium and technetium), 

stars and black holes). Together, these processes drive A , . , . „ , , „„, 

' , ° ' . AGB stars withm a certain mass range can thus be turned 

galactic chemical evolution and star formation. . , v , , . , 

° into carbon stars when the photosphcric carbon-to-oxygen 

Intermediate-mass stars, with Main-Sequence masses abundance ratio exceeds unity ( Be ssell, Wood & Lloyd 

of Mzams ~l-8 Mq, become Asymptotic Giant Branch g vans 1933) 

(AGB) stars before ending as a white dwarf on timescales . . . 

r ij.or.TVT 1 in/-iT- Massive stars also experience intense mass loss before 

Irom only t ~30 Myr up to as mucn as 10 Gyr. Luminous, 

r ,„4 T , , m ornn okaat^ j- they end through a supernova, leaving a neutron star 

Lagb ~ 10 4 L©, and cool, T c ff agb ~ 2500-3500 K, stars £ 6 ^ ' , 6 

, „ , n ' ■ • ; or black hole or deflagrating entirely on timescales 01 

on the upper slopes ol the AGB experience intense mass r\p 1 1 . . 

, n . „„ r ,, . . . , t ~1 to 30 Myr. Ol these the vast maiontv, stars with 

loss, losing up to 80 per cent ol their original mass at , r „ „_ , . , . , . 

rates of M AGB ~ lO^-lO" 4 M yr" 1 (Wood, Bessell & MzAMS 7 , \ P T S f supergiant (RSG 

Fox 1983; Wood et al. 1992). Their luminosity is normally f a § e ' ^ en they ^3™°^' ^ ~ 10 p 

, . , , , . , „ . , , , r Lp), and become cool, J c ff rsg ~ 3500-4000 K (Massey & 
generated m a hydrogen- burning shell, isolated from the » , ' . . . v 

, , j xi. 1 • Olscn 2003). As red supergiants they loose mass at rates 

degenerate carbon-oxygen core and the overlying convec- „ • ,„K,„an., i, 

6 J 6 of M - 10- J -10~ 3 Mq yr" 1 , but there is a lot of un- 

Send offprint requests to: e-mail: jaccoc3astro.keele.ac.uk certainty about the total amount of mass lost during the 

* based on observations obtained at the European Southern RSG phase, and therefore about the exact nature of the 

Observatory (La Silla, Chile) immediate supernova progenitors (Filippenko 1997). 
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AGB and RSG mass loss is essentially driven by the 
same mechanism: radiation pressure on dust grains which 
condense in the elevated atmospheres of radially pulsating 
giant stars and drag the gas along with it through mu- 
tual collisions. It seems therefore logical that the mass- 
loss rate would depend on the stellar luminosity (radia- 
tion pressure) and temperature (to allow dust to form). 
Indeed, both theory and observations confirm that this is 
the case. However, even the most sophisticated, hydrody- 
namical models suffer from a range of assumptions and 
are so far only presented for carbon stars, at an initial so- 
lar metallicity (e.g., Wachter et al. 2002; Sandin & Hofner 
2003). Mass-loss rates have been measured for a sample 
of AGB stars and red supergiants in the Large Magellanic 
Cloud (LMC) by van Loon et al. (1999), who demonstrate 
that the mass-loss rate increases with increasing luminos- 
ity, and for increasingly late spectral type for oxygen-rich 
(M-type) AGB stars and red supergiants, but they do not 
quantify these dependencies. 

We here present a recipe for the mass-loss rate as a 
function of luminosity and stellar effective temperature, 
M(L, T ff), for oxygen- rich AGB stars and red supergiants 
in the LMC, which undergo strong regular pulsations 
(Whitelock et al. 2003) and are enshrouded in dusty en- 
velopes. This is the first empirical relationship of its kind, 
that can be readily applied in models for stellar and galac- 
tic evolution. Although much work has been done on the 
mass-loss rates of cool giants in globular clusters and in the 
solar vicinity (e.g., Reimers 1975a, b; Jura & Kleinmann 
1992; Catelan 2000), these studies concentrate on first as- 
cent RGB stars and low-mass AGB stars with moderate 
mass-loss rates. Our formula is the first attempt to quan- 
tify the relationship between the measured values for the 
mass-loss rate, bolometric luminosity and stellar effective 
temperature for dusty cool giants. Valid well into the su- 
perwind regime of extremely high mass-loss rates, our for- 
mula provides a more direct alternative to the commonly 
used extra-polations of relationships which were derived 
for less massive and less evolved red giants. 

As part of this investigation, we obtained optical spec- 
tra for a number of highly obscured objects for which no 
such information was available, to estimate their stellar 
effective temperatures. We compiled infrared (IR) pho- 
tometry to sample the spectral energy distributions which 
were then modelled with a dust radiative transfer code 
in order to derive the mass-loss rates. We also present a 
more qualitative analysis of a small sample of carbon stars 
and (M)S-type stars with a carbon-to-oxygen ratio slightly 
less than unity, as well as a comparison with galactic AGB 
stars and red supergiants. 

2. Observational data 

2.1. Spectroscopy 

Spectra of (candidate) dust-enshrouded AGB stars and 
red supergiants were obtained on three occasions (1995, 
1998 and 2000). DFOSC at the 1.5m Danish telescope at 



Table 1. Spectroscopic targets, listed in order of increas- 
ing Right Ascension (J2000 coordinates). The last column 
refers to the observing run, where "D5" = DFOSC 1995, 
"E8" = EMMI 1998 and "E0" = EMMI 2000. 
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La Silla, Chile, was used from 8 to 11 December 1995 to 
obtain low-resolution (R ~ 900) long-slit spectra between 
~0.41 and 0.68 fim with grism #7 and a 1.5" slit width. 
EMMI at the ESO 3.5m NTT at La Silla, Chile, was used 
on 24 December 1998 and on 15 & 16 January 2000, to 
obtain low-resolution (R ~ 500) long-slit spectra between 
~0.61 and 1.04 /im with grism ^=4 and a 1.5" slit width. 
Integration times ranged between 5 and 30 minutes. 

The data were reduced using standard procedures in 
ESO-Midas. The frames were corrected for the electronic 
offset (bias) and the relative pixel response (Hatfield). 
Wavelength calibration was performed relative to He+Ar 
or He+Ne lamp spectra. The sky-subtracted spectra were 
corrected for the wavelength dependence of the instrumen- 
tal response as measured from the spectrum of a standard 
star, and for atmospheric continuum extinction. 

The targets (Table 1) were chosen from samples of 
mid-IR sources with known optically visible counterparts 
(dominating the DFOSC 1995 run) — mainly from Loup 
et al. (1997), bright mid-IR sources associated with heav- 
ily dust-enshrouded objects (mostly during the EMMI 
1998 run) — cf. Zijlstra et al. (1996), and red objects 
that were likely to be carbon stars (EMMI 2000 run). The 
latter were selected from Cioni et al. (2000) to add to 
the previous samples that turned out to be predominantly 
oxygen-rich. Several galactic M giants were also observed 
for reference (see Appendix A). 
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Table 2. Spectral classification and near-IR magnitudes (wavelengths in /zm). 
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2.2. Infrared photometry 

In order to estimate bolomctric luminosities and mass- 
loss rates, infrared photometry was compiled to sample 
the spectral energy distribution (SED). Near-IR photom- 
etry (Table 2) was retrieved from the DENIS catalogue 
(IJK S bands at 0.79, 1.23 and 2.15 ^m; Cioni et al. 2000; 
cf. Fouque et al. 2000) and the 2MASS catalogue (JHK S 
bands at 1.24, 1.66 and 2.16 ^tm; cf. Cohen, Wheaton 
& Megeath 2003). Mid-IR photometry (Table 3) was re- 
trieved from the MSX catalogue or from the original MSX 
images, from ISO (Trams et al. 1999), and from IRAS 
scans following the procedure described in Trams et al. 
(1999). Upper limits were determined where these would 



provide useful constraints on the shape of the SED, noting 
that the MSX upper limits correspond to 1-tr levels only. 

2.3. Additional samples 

We complement our sample with objects from van Loon 
et al. (1998) and Trams et al. (1999) for which spec- 
tral subclassification and mid-IR photometry is avail- 
able. No more carbon stars were included in this way, 
but several M-type and two MS-type stars were added. 
The MS star SHV 0524565-694559 was detected in the 
ISO mini-survey of the LMC (Loup et al., in prepara- 
tion), with F 4 . 5 = 16.2 ± 2.3 mJy and F 12 = 5.3 ± 0.7 
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Table 3. Mid-IR flux densities (in Jy, with the wavelengths in /im). The column "s" indicates the availability of an 
ISO spectrum ("y"). Values accompanied by a colon are uncertain, whilst question marks indicate detections that are 
suspect. 
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mJy. Another possible S star listed in Trams et al. is 
SHV 0522118— 702517, but its classification is uncertain 
and no subclassification is given. 



A further 34 IRAS-detected objects with spectral sub- 
classification are available from Loup et al. (1997; their 
Table 1); these are predominantly early-M type super- 
giants that are already well represented in our sample. 



3. Analysis 

3.1. Spectral classification 

The spectral type of oxygen-rich (M-type) stars was de- 
rived from the relative strengths of the TiO and VO molec- 
ular absorption bands, and of the Ca II triplet, whilst the 
spectral type of carbon-rich (C-type) stars was derived 
from the relative strengths of the C2 and CN molecular ab- 
sorption bands. S-type stars with a carbon abundance ap- 
proaching that of oxygen have prominent ZrO bands, with 
MS-type stars representing a milder form of this chemical 
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Fig. 1. DFOSC spectra of M-type supergiants. 



peculiarity. For more details about these bands the reader 
is referred to Turnshek et al. (1985). 
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Fig. 2. EMMI spectra of M-type AGB stars, displayed 
twice on different scales to exploit the full dynamic range 
of the data. 



3.1.1. M-type stars 

The M-type stars observed with DFOSC (Fig. 1) turn out 
to be red supergiants, whilst the M-type stars observed 
with EMMI (Fig. 2) are on the AGB (see Table 4, Fig. 
11). The 4000 to 7000 A region is not optimal for classifi- 
cation of late-M types. The 7000 to 9000 A region is much 
better, chiefly due to the diagnostic value of the VO bands. 
There can be some ambiguity, notably in WOH G64 (Fig. 
1) where various band strengths (in particular at 6200 



A) suggest a very late type (M7-8) whilst some relative 
band strengths (e.g., 6200 compared to 5900 A) suggest 
a much earlier type (~M5). Band strength ratios might 
better reflect the excitation temperature, whilst individ- 
ual bands may appear stronger as a result of enhanced 
optical depth in the puffed up atmospheres of mass-losing 
giants. Nonetheless, normally we obtain a spectral classi- 
fication that is accurate enough for our purposes. 
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Fig. 4. EMMI spectra of carbon stars. 
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Fig. 3. DFOSC 
GRV 0519-6700. 



spectrum of the S-type star 



3.1.2. An S-type star 

In S-type stars, the 5469 A VO band becomes weaker (e.g., 
relative to the ZrO bands) with increasing C/O abundance 
ratio. It is quite strong in GRV 0519-6700 (Fig. 3), which 
must therefore have a C/O ratio not too close to unity. The 
YO bands at 5972 and 6132 A, and the pair of roughly 
symmetrical CaCl bands around 6200 A only appear in the 
cooler S-type stars (Keenan & Boeshaar 1980), and they 
are clearly present in the spectrum of GRV 0519— 6700. 
There is no sign of 13 CN bands suppressing the peak in the 
spectrum between 6450 and 6456 A, and one can assume 
that this star is not particularly 13 C rich. 
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3.1.3. Carbon stars 

The spectra of carbon stars (Figs. 4 & 5) are generally 
quite similar, with a few exceptions. The C/O abundance 
ratio determines the relative strength of the C2 and CN 
bands. The C2 bands are best measured shortward of 5600 
A, but unfortunately that spectral region was not included 
(and the stars are much fainter in the V band). The best 
diagnostic in our spectral coverage is the collection of C2 
bands shortward of 6200 A, and around 7700 and 8800 A. 
The temperature is best measured from the strength of 
the CN bands, roug hly between 6300 and 7500 A. 

DCMC 050738.74-733233.4 (Fig. 4) seems to be a 
warm carbon star and may either be extrinsic (i.e., en- 
riched by material from a previous AGB companion), ex- 
periencing the effects of a thermal pulse, or on a post- AGB 
track. As a consequence of the relatively high temperature 
all bands become weak and it is therefore impossible to es- 
timate the C/O ratio from our spectrum. 

IRAS 05289-6617 (Fig. 5) exhibits extremely strong 
absorption in some of the CN bands, with no sign of C2. 
This dust-enshrouded object might be nitrogen-enriched, 
possibly as a result of (mild) HBB. The mid-IR ISOCAM- 
CVF spectrum (Trams et al. 1999) shows an emission 
feature around 11.5 (im that could either be due to SiC 
grains or Polycyclic Aromatic Hydrocarbons (PAHs). This 
star is located at the fringes of the populous elliptical 
intermediate-age cluster NGC 1978 and its colours suggest 
the lack of warm dust (van Loon et al., in preparation). 
The latter might be explained if the star has just recovered 
from a thermal pulse. 

3.1.4. Emission-line objects 

Two dusty objects turned out to exhibit a rich spectrum 
of strong emission lines (Fig. 6) including Balmer Ha, the 
Ca II triplet and the Paschen hydrogen series longward of 
8400 A. There is evidence for a two-component spectrum, 
with a warm (or hot) continuum peaking around 6000 A 
and a cold continuum dominating longward of 9000 A. The 
cool component seems to be relatively more important in 
IRAS 04530-6916 than in LI-LMC 1522. 

IRAS 04530-6916 has strong [S 11] emission in the 
6717+6731 A doublet. This very luminous object was in- 
terpreted as a candidate massive Young Stellar Object on 
the basis of its IR colours (van Loon et al. 2001b). 

Besides the broad Ha line, LI-LMC 1522 shows strong 
He 1 emission lines at 6678 & 7065 A as well as emission 
in [Ca 11] at 7291+7324 A. Strong He lines are also seen 
in objects such as ryCarinae (Wallerstein et al. 2001), and 
LI-LMC 1522 may be an evolved supergiant. Indeed, its 
spectrum resembles that of the B[e] supergiant HenS134 
and like that luminous B0 star also shows the 6170 A TiO 
band in emission and numerous Fe 11 emission lines in the 
6000-6500 A region (Lamcrs et al. 1998). 

When taking the spectrum of the dust-enshrouded 
AGB star IRAS 05003-6712, the B2 emission-line star 
Al 50 (Andrews & Lindsay 1964) happened to fall on the 
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Fig. 5. EMMI spectrum of the carbon star 
IRAS 05289-6617. 
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Fig. 7. DFOSC spectrum of the emission-line object Al 50. 



slit of the spectrograph too. Its spectrum is characterised 
by bright emission lines from the hydrogen Balmer series 
and helium, superimposed on a hot continuum (Fig. 7). 
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Table 4. Results from the modelling of the spectral en- 
ergy distributions: the stellar effective temperature (from 
the optical spectrum), T e g, the temperature at the inner 
radius of the dust envelope, Td ust , the wind speed, v, the 
bolometric luminosity, L, and the total (gas+dust) mass- 
loss rate, M. We omitted the DCMC prefix to the first 
four carbon stars. 
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3.2. Modelling the spectral energy distributions 

The spectral energy distribution (SED) of each object 
was reproduced with the dust radiative transfer model 
DUSTY (Ivezic, Nenkova & Elitzur 1999), which was made 
to calculate the radial density distribution consistent with 
radiation-driven wind theory (Ivezic & Elitzur 1995). This 
takes into account the dependence of the wind speed on 
the stellar parameters luminosity and metallicity (via the 
dust-to-gas mass ratio), and the drift speed between the 
dust and the gas. This formalism is confirmed to hold 



for OH/IR stars in the LMC for which the wind speeds 
have been measured from the OH maser emission profiles 
(Marshall et al. 2004), a few of which are in common with 
the sample we discuss here. Some of the model parameters 
that we will analyse in detail are summarised in Table 4. 
We discuss some of the input parameters and the reliabil- 
ity of the model output, before analysing the results. 



3.2.1. Stellar effective temperatures 

The stellar effective temperatures corresponding to the M- 
type spectral subclasses are taken from Fluks et al. (1994). 
These are strictly speaking valid for luminosity class III 
giants of solar abundances, and differences may be ex- 
pected for metal-poor stars and for luminous supergiants. 
The stratification of the extended molecular atmosphere 
together with variations over the pulsation cycle limit the 
accuracy with which these cool giants may be assigned 
a value for the effective temperature. Where more de- 
terminations of the spectral type were available and/or 
classifications to half a subclass were given, the spectral 
subclasses were averaged and the stellar effective temper- 
ature was obtained by linear interpolation of the given 
values. The thus adopted values are probably only accu- 
rate to ~300 K (cf. Houdashelt et al. 2000). We used the 
synthetic spectra from Fluks et al. closest to the actual 
object's spectral type as models for the external radiation 
field when running the DUSTY code. 

Despite the prominent and unique molecular absorp- 
tion features characterising carbon star spectra, it is re- 
markably difficult to classify these spectra and even more 
difficult to assign stellar effective temperatures. We con- 
structed a linear grid of stellar effective temperature val- 
ues between 2200 K for spectral type C10 and 4000 K 
for CO. We stress that this is a purely ad hoc tempera- 
ture scale, and we will limit the quantitative interpreta- 
tion of the carbon stars in our sample, but it nevertheless 
allows to rank the carbon stars according to their tem- 
peratures and study the qualitative effect that differences 
in the surface temperature may have on the circumstcllar 
envelope and mass- loss rate. Bergeat, Knapik & Rutily 
(2001) construct a photometric temperature classification 
scheme, and it would be useful to tally their findings with 
a spectral classification scheme. In running dusty we as- 
sumed a blackbody of the effective temperature as the 
model for the external radiation field. 

For MS and S-type stars things become even more 
uncertain. We assign stellar effective temperatures in be- 
tween the values we would have given were they of either 
M-type or C-type. For the purpose of running DUSTY, we 
assumed that the spectral energy distribution of the ex- 
ternal radiation field were that of an equivalent M-type 
synthetic spectrum from Fluks et al. (1994). 
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3.2.2. Grain properties 

Because of the limitations of the data to constrain the dust 
grain properties (size, shape, composition...), we have re- 
stricted ourselves in the choice of dust type, and assumed 
a grain size of 0.1 /im unless the fit of the dusty model to 
the observed SED improved notably by making a different 
choice. In order to derive the total mass budget one also 
has to adopt a mass density for the grains. Although this 
varies with the type of material the grains are composed 
of, we adopt a default value of /9 gra in = 3 g cm~ 3 . 

For oxygen-rich M-type stars and MS and S- 
type stars we adopted the astronomical silicate from 
Draine & Lee (1984), but for the carbon stars 
we took different mixtures of graphite (Draine & 
Lee 1984), amorphous carbon (Hcnning & Mutschke 
1997) and silicon carbide (SiC, from Pegourie 1988). 
By default we used amorphous carbon only, but 
DCMCJ050738.74- 733233.4, SHV 0500533-690617 and 
SHV 0504353-712622 required the inclusion of 10 per 
cent of SiC, and the ISO spectrum of IRAS 05289-6617 
suggests the need of 20 percent of SiC (van Loon ct 
al. 1999). The SEDs of DCMC J050024.95-710130.7 and 
DCMCJ050615. 14- 720657.1 are better reproduced by 
taking only 50 per cent of amorphous carbon, with the 
remainder made up with graphite or with 40 per cent 
graphite and 10 per cent SiC, respectively. 

The near-IR extinction and mid-IR emission of several 
amongst the most dust-enshrouded M-type stars of our 
sample could be better brought in line with each other 
by reducing the grain size to 0.06 ^m (IRAS 5128—6455) 
or 0.05 /xm (IRAS 04407-7000, IRAS 04498-6842 and 
WOH G64) , or by taking a standard MRN size distribu- 
tion (a power-law with exponent —3.5; Mathis, Rumpl & 
Nordsieck 1977) between a minimum grain size of 0.01 /im 
and a maximum grain size of 0.1 /xm (IRAS 05003—6712). 
Amongst the carbon stars we deviated from the 0.1 ^m 
grain size only in the case of SHV 0504353-712622 (0.01 
H, and IRAS 5289-6617 for which we adopted an MRN 
distribution of larger grains between 0.1 and 1 /j,m. The 
latter object has both a cool photosphere and abnormally 
cool dust, and with the high fraction of SiC (see above) it 
is also mineralogically peculiar. 

Models were first computed for a dust temperature at 
the inner radius of the circumstellar envelope of Td us t = 
1000 to 1200 K, but in many cases an acceptable fit to 
the observed SED could only be achieved by lowering this 
temperature significantly — although altering the grain 
properties may also have the desired effect. 

3.2.3. Luminosities 

To derive the bolometric luminosities the DUSTY output 
SED needs to be scaled to fit the observed SED. Adopting 
a distance to the LMC of 50 kpc, this scaling factor then 
yields the bolometric luminosity under the assumption 
that interstellar extinction is negligible over most of the 
infrared SED. 



All objects are long-period variables. Although the am- 
plitude of variability drastically diminishes at IR wave- 
lengths, it can still reach a magnitude in bolometric lu- 
minosity for AGB stars (for supergiants the amplitude 
is usually smaller in units of magnitude even if it cor- 
responds to a larger change in energy output). This vari- 
ability is reflected in differences between the DENIS and 
2MASS near-IR magnitudes, and between the MSX, ISO 
and IRAS mid-IR flux densities. We attempt to fit the 
dusty model to a roughly "mean" SED. Apart from the 
uncertainty in the distance to the LMC of around 10 per 
cent, the values for the luminosities of individual stars are 
believed to be accurate to within ^20 per cent. 



3.2.4. Mass-loss rates 

In order to obtain total (gas+dust) mass-loss rates one 
has to know the gas-to-dust mass ratio in the outflow. 
In the absence of knowledge thereof we assume a value 
of ip = Pgas/Pdust = 500, which we justify by scaling the 
typical value of ip = 200 for giants with solar abundances 
to the lower metallicity of the LMC (cf. van Loon et al. 
1999). For most of the dust-enshrouded objects it would 
seem plausible that the dust condensation process reaches 
its maximum efficiency, supporting these low gas-to-dust 
ratios, but for less evolved and/or hotter objects this may 
no longer be the case, possibly leading to underestimated 
values for their mass-loss rates. 

Many of the objects in our sample are detected in the 
IRAS data at 60 /im, at a level much higher than can be 
reconciled with a constant mass-loss rate model. Although 
these detections need to be confirmed at a higher sensi- 
tivity and angular resolution, possible explanations could 
include the heating of interstellar dust by the supergiants 
and/or ancient episodes of enhanced mass loss. The lat- 
ter is expected especially for AGB stars descending from 
low-mass (1-3 M ) Main-Sequence progenitors, as the ef- 
fects of thermal pulses on the upper AGB on the surface 
temperature and luminosity are much more dramatic than 
for their more massive siblings (cf. Schroder, Winters & 
Sedlmayr 1999). 

There is ample opportunity for making different 
choices of grain properties, dust temperature, optical 
depth etcetera without affecting the fit of the DUSTY 
model to the SED. This is especially true for objects with 
optically thin dust shells. However, fortunately the result- 
ing mass-loss rate is surprisingly little influenced by this 
degeneracy, thanks to the scaling inter-relationships be- 
tween the quantities that govern the radiative transfer 
(Elitzur & Ivezic 2001; cf. van Loon et al. 1997). Apart 
from largely systematic uncertainties such as the gas-to- 
dust ratio, the values for the mass-loss rates of individual 
objects are reliable to within a factor of ~two. This is con- 
firmed by the observed spread in the L, T c g dependence 
of the mass-loss rate (see below). 



10 



van Loon, Cioni, Zijlstra, Loup: Mass-loss rates of red supergiants and AGB stars 



o ^-t) 



8 o 



o o 



JO o 



AGB RSG 



O = M stars 
a = MS/S star; 
• = C stars 

J III L 



3.5 



4.5 5 5.5 

log L [L s ] 



Fig. 8. Mass-loss rate as a function of bolometric lumi- 
nosity, for M-type stars (circles), MS or S-type stars (tri- 
angles), and carbon stars (solid dots). The classical AGB 
luminosity limit is marked by a vertical dotted line. The 
long and short-dashed slanted lines mark the classical and 
multiple-scattering limits to the mass-loss rate (van Loon 
et al. 1999). 

4. Mass-loss rate as a function of bolometric 
luminosity and stellar effective temperature 

4.1. Mass-loss rates from oxygen-rich dust-enshrouded 
red giants in the Large Magellanic Cloud 

The mass-loss rates reached by more luminous AGB stars 
are higher than those reached by less luminous AGB stars, 
but there is a lot of intrinsic scatter in the mass-loss rate 
versus luminosity relation due to evolutionary effects (Fig. 
8). Stars that have only just reached the thermal pulsing 
stage of AGB evolution do not experience as high a mass- 
loss rate as when they reach the AGB tip and develop 
a superwind. Likewise, although red supergiants feature 
mass-loss rates that are comparable to those reached by 
luminous tip- AGB stars, they do not seem to continue the 
trend of the mass-loss rate increasing with luminosity be- 
cause their true superwind stage lasts very briefly and thus 
this evolutionary phase is largely missed. The one excep- 
tion within our sample is WOHG64, the most luminous 
red supergiant in the Magellanic Clouds with the largest 
mass-loss rate. 

Our sample is obviously biased against the most 
severely dust-enshrouded stars for which no spectrum of 
the stellar photosphere can be taken at optical wave- 
lengths. This is especially true for the carbon stars, of 
which the circumstellar dust envelopes become optically 
thick relatively easily as a result of the higher opacity of 
carbon-rich dust and the more compact geometry of the 
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Fig. 9. Mass-loss rate as a function of stellar effective tem- 
perature. Symbols are the same as in Fig. 8. 



envelope (cf. Fig. 13). Hence our sample seems to have 
fewer examples of the extremely high mass-loss rates as 
compared with the sample of obscured objects in van Loon 
et al. (1999) and their relationships for the classical mass- 
loss rate limit and mass-loss rate limit in the presence of 
multiple scattering of photons off dust grains (the long and 
short-dashed lines in Fig. 8, respectively). Nonetheless, 
our sample does include several OH/IR stars which are 
in the superwind phase, and our optical spectroscopy of 
these heavily dust-enshrouded objects makes a significant 
contribution to probing the stellar parameters well into 
this literally obscure phase of stellar evolution. 

At least part of the evolutionary effect in the disper- 
sion of mass-loss rates can be related to the stellar photo- 
spheric temperature. On the AGB, stars evolve not only to 
become brighter but also to become cooler, whilst red su- 
pergiants first cool rapidly but are then expected (on the 
basis of models) to increase in luminosity somewhat at an 
essentially constant temperature. It is not clear at what 
precise evolutionary stage we capture the red supergiants 
in our sample. 

When the mass-loss rate is held against the stellar ef- 
fective temperature, however, no clear trend is seen (Fig. 
9). The M-type stars with the lowest mass- loss rates are 
fairly warm and the coolest M-type stars have relatively 
high mass-loss rates, but the red supergiants trouble this 
picture as they have high mass-loss rates not because they 
are cool (they are, in fact, amongst the warmer stars in 
our sample) but simply because they are luminous. This 
might give the impression that the luminosity dependence 
of the mass-loss rate dominates over the temperature de- 
pendence. However, it should be realised that the range 
in luminosities is vastly larger (two orders of magnitude) 
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than the range in stellar effective temperature (less than 
a factor two) within our sample. 

A surface can be fit to the M(L, T e g ) plane for the M- 
type stars, with useful constraints on the fit parameters. In 
practice, we fit a flat surface in logarithmic space using a 
non-linear least-squares Marquardt-Levenberg algorithm. 
We thus derive: 

log M = a + ^{ ids^) +7 log GSk)' 

with for M-type stars: a = -5.65 ± 0.15, /3 = 1.05 ± 0.14 
and 7 = -6.3 ± 1.2. 

This solution for the mass-loss rate recipe features 
an approximately linear increase of the mass-loss rate 
with an increase in luminosity, which is consistent with 
the radiation-driven dust wind theory including multiple- 
scattering (cf. van Loon et al. 1999). This is also repro- 
duced for a set of synthetic stars by Schroder, Wachter & 
Winters (2003), who find a luminosity dependence with 
an exponent of 1.16. 

Besides the luminosity dependence, the mass-loss rate 
is also seen to depend sensitively on the stellar effective 
temperature, approximately inversely as the sixth power. 
A strong temperature dependence was initially suggested 
by Arndt, Fleischer & Sedlmayr (1997) on the basis of 
models for carbon stars. Although they found an even 
stronger dependence of the mass-loss rate on the stellar 
effective temperature, with an exponent of —8.3, their up- 
dated model calculations reduce the temperature expo- 
nent to —6.8 (Wachter et al. 2002); this is fully consistent 
with our observed value of —6.3 ± 1.2. 

These authors also found a luminosity dependence of 
the mass- loss rate, with an exponent of 1.5 (Arndt et al. 
1997) and 2.5 (Wachter et al. 2002), respectively. Although 
these values are significantly larger than our observed 
value for the luminosity exponent of 1.05 ± 0.14 they also 
include a dependence on the stellar mass, with a decreas- 
ing mass-loss rate for an increase in stellar mass. This 
effect counteracts upon the luminosity dependence and, if 
the mass dependence were somehow incorporated within 
the luminosity dependence, it would support a lower ob- 
served value for the luminosity exponent in the absence 
of an explicit parameterisation in terms of stellar mass. 
Our formula may thus be consistent with the models of 
Arndt and Wachter, but this must be tested with objects 
for which the current masses are determined from accurate 
modelling of their pulsation characteristics. 

The residuals from the M(L, T ? s) fit (Fig. 10) show a 
dispersion of only 0.3 dex in logAf (a factor two in M), 
which is comparable to the accuracy of an individual esti- 
mate of the mass-loss rate. Other stellar parameters may 
contribute to the spread, but it is clear that the luminos- 
ity and temperature dominate the mass-loss rate recipe 
for these stars. 

We do not account for departures from spherical sym- 
metry of the mass loss. Asymmetric flows are often seen 
in Planetary Nebulae and post-RSG objects, and thought 
to have their origin during the previous superwind phase. 
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Fig. 10. Discrepancies between measured mass-loss rates 
from modelling of the SED, and the mass- loss rates as 
predicted from the L, T c g dependence of the mass-loss 
rate for M-type stars, plotted versus bolometric luminosity 
(a) and stellar effective temperature (b) . Symbols are the 
same as in Fig. 8. 



Much amplified in post-AGB and post-RSG objects due 
to the interaction of a fast wind with the slower super- 
wind, it is not yet clear how prominent or common these 
asymmetries are in the winds of AGB stars and red super- 
giants. There exists no evidence for the global mass-loss 
rate to be affected by asymmetry in the outflow kinemat- 
ics, or by the presence of a binary companion. Our sample 
might already include binaries or objects with asymmet- 
ric outflows as we have not wittingly excluded them — in 
which case the mass-loss rates for such objects cannot be 
different by more than a factor two. 

We now divide the sample of M-type stars into red su- 
pergiants, with logi > 4.9, and AGB stars, with logL < 
4.9 (see next section), and fit a similar M(L,T c g) for- 
mula to each of these subsamples. Despite the smaller 
number of stars to constrain the fit (10 red supergiants 
and 12 AGB stars), the results are in excellent agreement 
with each other. In particular, the luminosity exponent 
/3 = 0.82±0.38 for the red supergiants and (3 = 1.10 ±0.40 
for the AGB stars are indistinguishable. Despite the much 
larger uncertainty on the temperature exponent, where 
we note that the distribution of red supergiants is heavily 
skewed towards early types, with J — — 10.8±2.7for the 
red supergiants and 7 = —5.2 ± 3.0 for the AGB stars 
these are also fully consistent with each other. The fits 
lead to an estimate for the reference mass-loss rate of 
a = —5.3 ± 0.5 and a = —5.6 ± 0.2 for the red super- 
giants and AGB stars, respectively, again suggesting that 
the mass-loss rates from both red supergiants and oxygen- 
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rich AGB stars are governed by the same mechanism and 
in the same manner. 

If the mass-loss rate recipe for M-type stars is applied 
to the (M)S and C-type stars, their dispersion is reduced 
as well, suggesting that the mass-loss rate of these types 
of chemically distinct objects may depend on the luminos- 
ity and stellar effective temperature in a similar fashion. 
Attempts to fit an M(L,T e s) relation to these classes of 
stars separately lead to qualitative similar exponents (pos- 
itive for the luminosity dependence and negative for the 
temperature dependence) but in view of the smaller num- 
ber of stars and the uncertainty with regard to the tem- 
perature scale for these objects we prefer not to publish 
the results of these fits to avoid abuse. 

4.2. Comparison with galactic red giants 

We now investigate to what extent the mass-loss rate 
recipe is valid for other samples, by applying it to galac- 
tic AGB stars and red supergiants. Jura & Kleinmann 
published studies of the mass loss from red giants in the 
solar neighbourhood, which provides us with three sam- 
ples: (1) AGB stars pulsating with an intermediate-period 
(Jura & Kleinmann 1992), (2) dust-enshrouded AGB stars 
with a large 60 /im excess (Jura & Kleinmann 1989), and 
(3) red supergiants (Jura & Kleinmann 1990). From their 
spatial distributions Jura & Kleinmann estimate progeni- 
tor masses of respectively 1-1.2 M Q , 1.5 M Q and >20M Q 
for these three samples. We selected all M-type stars for 
which spectral sub-classification is available to allow a de- 
termination of T e ff, and only those stars for which Jura 
& Kleinmann estimated a mass-loss rate: 5 intermediate- 
period AGB stars (with pulsation periods of 3 1 0-372 d), 
26 dust-enshrouded AGB stars and 19 red supergiants. 
Their properties are listed in Table 5. Jura & Kleinmann 
(1990) list spectral types for the supergiants, but spectral 
types were obtained from Simbad for the other two sam- 
ples and for the red supergiant NML Cyg. The values for 
T e ff are assigned according to the same calibration scheme 
as for our sample of magellanic objects. 

Jura & Kleinmann (1992) estimated distances for the 
intermediate-period AGB stars using the K-band period- 
luminosity relation for Mira variables. For a few nearby 
examples these are consistent with determinations from 
Hipparcos data (Knapp et al. 2003). Nevertheless, Jura & 
Kleinmann assumed the luminosities to be 5,000 Lq. They 
estimated distances for the dust-enshrouded AGB stars 
from the bolometric fluxes, where they assumed that they 
all have a luminosity of 10 4 L Q (Jura & Kleinmann 1989). 
For the supergiants they estimated distances in most part 
from cluster membership (Jura & Kleinmann 1990). The 
distance ranges for the samples are 80-300 pc, 270-1 100 pc 
and 200-2500 pc, respectively. 

The method of Jura & Kleinmann to derive mass-loss 
rates, which they applied to all three samples in the same 
manner, is based on a scaling of the measured 60 /jm flux 
density (Jura 1986, 1987). The result depends on distance, 



Table 5. Galactic samples of red giants (Jura & 
Kleinmann 1989, 1990, 1992), with the spectral type, the 
stellar effective temperature, T e g, the bolometric luminos- 
ity, L, and the mass- loss rates as determined by Jura & 
Kleinmann, M 1; and as predicted by our recipe, M 2 . We 
omitted the IRAS prefix to some of the dust-enshrouded 
AGB stars. 



Object Spectral T eS logL log Mi logM 2 





Type 


K 


Lq 


M Q /yr 


M /yr 


Intermediate-pe 


riod AGB stars 


(Jura & Kleinmann 1992) 


SCrB 


M7 


3126 


3.70 


-6.96 


-5.66 


RLeo 


M8 


2890 


3.70 


-7.07 


-5.44 


RLMi 


M7 


3126 


3.70 


-6.82 


-5.66 


RXBoo 


M7.5 


3008 


3.70 


-6.80 


-5.55 


WHya 


M7 


3126 


3.70 


-7.00 


-5.66 


Dust- enshrouded AGB stars (Jura & Kleinmann 1989) 


KUAnd 


M9 


2667 


4.00 


-4.72 


-4.91 


WXPsc 


M9 


2667 


4.00 


-4.64 


-4.91 


02316+6455 


M9 


2667 


4.00 


-5.07 


-4.91 


02351-2711 


M9 


2667 


4.00 


-5.44 


-4.91 


NML Tau 


M6 


3309 


4.00 


-5.22 


-5.50 


04307+6210 


M6 


3309 


4.00 


-5.44 


-5.50 


TXCam 


M8.5 


2778 


4.00 


-5.60 


-5.02 


05411+6957 


M9 


2667 


4.00 


-5.23 


-4.91 


VCam 


M7 


3126 


4.00 


-5.66 


-5.34 


06300+6058 


M9 


2667 


4.00 


-5.28 


-4.91 


GXMon 


M9 


2667 


4.00 


-5.03 


-4.91 


IWHya 


M9 


2667 


4.00 


-5.02 


-4.91 


V2108 Oph 


M9 


2667 


4.00 


-5.42 


-4.91 


V774 Sgr 


M5 


3434 


4.00 


-5.55 


-5.60 


18009-2019 


M8 


2890 


4.00 


-5.37 


-5.12 


18135-1641 


M5 


3434 


4.00 


-5.30 


-5.60 


18204-1344 


M8 


2890 


4.00 


-5.38 


-5.12 


Villi Oph 


M9 


2667 


4.00 


-5.28 


-4.91 


18413+1354 


M7 


3126 


4.00 


-5.26 


-5.34 


V3953 Sgr 


M9 


2667 


4.00 


-5.17 


-4.91 


V3880 Sgr 


M8 


2890 


4.00 


-5.00 


-5.12 


V342 Sgr 


M9 


2667 


4.00 


-5.36 


-4.91 


20440-0105 


M9 


2667 


4.00 


-5.57 


-4.91 


UUPeg 


M7 


3126 


4.00 


-5.80 


-5.34 


RTCep 


M6 


3309 


4.00 


-5.60 


-5.50 


23496+6131 


M9 


2667 


4.00 


-5.27 


-4.91 


Red supergiants 


(Jura & Kleinmann 1990) 




KWSgr 


MOIa 


3895 


5.00 


-5.40 


-4.89 


VXSgr 


M4Ia 


3574 


5.18 


-4.40 


-4.47 


UYSct 


M4Iab 


3574 


5.00 


-5.15 


-4.66 


BD+24°3902 


Mlla 


3810 


5.00 


-6.00 


-4.83 


BICyg 


M4Ia 


3574 


5.30 


-5.22 


-4.34 


KYCyg 


M3.5Ia 


3620 


5.60 


-5.40 


-4.06 


NML Cyg 


M6III 


3309 


5.70 


-4.00 


-3.71 


H Cep 


M2Ia 


3736 


5.60 


-6.05 


-4.15 


PZCas 


M2Ia 


3736 


5.30 


-5.00 


-4.46 


TZ Cas 


M2Iab 


3736 


5.00 


-5.52 


-4.78 


SPer 


M3Ia 


3666 


5.30 


-5.15 


-4.41 


SUPer 


M3Iab 


3666 


5.00 


-5.22 


-4.73 


a Ori 


Ml 


3810 


5.00 


-6.70 


-4.83 


VYCMa 


M5I 


3434 


5.60 


-4.00 


-3.92 


EVCar 


M4Ia 


3574 


5.30 


-5.22 


-4.34 


CKCar 


M3.5Iab 


3620 


5.00 


-5.70 


-4.69 


IX Car 


M2Iab 


3736 


5.00 


-5.05 


-4.78 


V396 Cen 


M4Ia 


3574 


5.30 


-6.00 


-4.34 


AHSco 


M4 


3574 


5.48 


-5.10 


-4.16 
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Fig. 11. Discrepancies between the mass- loss rates de- 
rived from 60 /im flux densities, and the mass- loss rates as 
predicted from the M(L, T e s) formula for dust-enshrouded 
stars in the LMC, plotted versus stellar effective tem- 
perature, for stars in the solar neighbourhood (Jura & 
Kleinmann 1989, 1990, 1992). 

d, luminosity, L, and wind speed, v, as M oc vd 2 /y/L. The 
wind speed was measured for almost every intermediate- 
period and dust-enshrouded AGB star that we selected, 
but for only 7 out of the 19 supergiants. For the stars with- 
out a direct measurement of the wind speed they assumed 
v = 15 and 30kms _1 for dust-enshrouded AGB stars 
and supergiants, respectively. The selected intermediate- 
period AGB stars have a (measured) wind speed of only 
v ~ 6kms~ 1 . Some supergiants also have a slow wind, 
notably aOri and fiCcp (v = lOkms -1 ) and Jura & 
Kleinmann (1990) suggest that this may be due to a much 
reduced dust-to-gas ratio (ip ^> 200). If that is the case 
then their values for the mass-loss rates of these stars will 
have been underestimated. 

There is excellent agreement between the mass-loss 
rates as predicted using our recipe and the measured 
Jura & Kleinmann mass-loss rates for the galactic dust- 
enshrouded AGB stars (Fig. 11, solid squares), where most 
of the scatter is due to deviations from the assumed lumi- 
nosity of 10 4 L for these stars. We note that the sample of 
dust-enshrouded AGB stars in Jura & Kleinmann (1989) 
includes another 5 more extreme oxygen-rich objects with- 
out spectral sub-classification, for which they estimate 
mass-loss rates between 2 x 10 -5 and 2 x 10~ 4 M Q yr _1 . 
A similar bias is present in our magcllanic sample (see 
next section). The good agreement between the dust- 
enshrouded (oxygen-rich) AGB stars in the LMC and in 
the solar neighbourhood suggests that there is no evidence 
for the total mass-loss rates for these classes of stars to de- 



pend on the initial metal abundance, a conclusion which 
was already reached by van Loon (2000). 

Our recipe also correctly predicts the mass-loss rates 
for the most extreme red supergiants, notably NML Cyg, 
VX Sgr and VY CMa that have mass-loss rates of the or- 
der 10 _4 M©yr _1 (Fig. 11, circles). However, most of the 
red supergiants from the Jura & Kleinmann (1990) sample 
have measured mass-loss rates that are significantly lower 
than the predicted values by about a factor 3-50 (although 
most are within a factor 10 of the predictions). This can 
be due to invalidity of the mass-loss recipe or the assump- 
tions made in the measured values. We already mentioned 
that the most discrepant objects in our comparison, a Ori 
and \i Cep have anomalously slow winds and probably a 
low dust-to-gas ratio. Furthermore, these two stars are the 
nearest supergiants in the sample and they are spatially 
resolved at 60 /im, which may have led to a further under- 
estimate of the mass-loss rate. The same may be true to 
a lesser extent for the other supergiants. 

There is a clear discrepancy for the intermediate- 
period AGB stars (Fig. 11, triangles), with the predictions 
exceeding the measured values by a factor 10-40. These 
are nearby objects and the 60 fim flux may again have been 
underestimated. They are cool and pulsate vigorously: 4 
out of 5 are Mira-type variables whilst WHya is a semi- 
regular variable. But their periods are with P < 400 d 
shorter than those of the dust-enshrouded stars, super- 
giants and the magcllanic sample. They are also optically 
bright. Our recipe is not meant to be applied to these 
stars as it was derived for dust-enshrouded objects with 
P > 400 d. The dust formation in the intermediate-period 
AGB stars may not be at its maximum efficiency, resulting 
in lower dust-to-gas ratios (explaining their slow winds). 
Perhaps the different pulsation properties or the higher 
gravity near the photospheres of these stars causes the 
atmosphere to be elevated less than in more extreme ob- 
jects. Indeed, circumstellar dust envelopes are much more 
prominent for objects with P > 400 d than for objects 
with P < 400 d (e.g., Knapp et al. 2003). 

5. The Hertzsprung-Russell diagram of 

dust-enshrouded AGB stars and red supergiants 

Now that we have determined the bolometric luminosi- 
ties and stellar effective temperatures for the objects in 
our sample, we can place them in a physical Hertzsprung- 
Russell diagram (Fig. 12). This allows a comparison with 
theoretical evolutionary tracks to be made. To this pur- 
pose, the tracks from Bertelli et al. (1994) are plotted, for 
a typical LMC metallicity of [Fe/H]=-0.4. 

Many of the stars in our sample are cooler than the 
endpoints on the AGB or RSG tracks. Indeed, despite the 
low metal abundance in the LMC, very late spectral types 
are encountered, with now at least two M10 stars known 
on the AGB and for instance the red supergiant WOH G64 
also being much cooler than the evolutionary models pro- 
duce (cf. Massey & Olsen 2003). The late-type AGB stars 
in our sample are not even the most extreme objects; more 
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Fig. 12. The Hcrtzsprung- Russell diagram of bolometric 
luminosity versus stellar effective temperature (symbols 
are the same as in Fig. 8). The classical AGB luminosity 
limit is marked by a horizontal dotted line. Three cool 
stars just above this limit may be AGB stars whose lu- 
minosity is enhanced as a result of Hot Bottom Burning 
(HBB). Evolutionary tracks from Bertelli et al. (1994) are 
plotted and labelled by their Main Sequence progenitor 
mass. 



heavily obscured OH/IR stars such as IRAS 05298—6957 
(Wood et al. 1992; van Loon et al. 2001a) exist, for which 
no optical spectrum has yet been possible to take. 

On the basis of the evolutionary tracks, the classical 
AGB luminosity limit and the distribution of the stars 
in the diagram, it seems a-posteriori fairly obvious which 
stars are intermediate-mass AGB stars and which are 
more massive supergiants. Rather than interpreting the 
three cool stars just above the classical AGB limit (Fig. 
12) as much cooler than expected supergiants, it seems 
more logical to interpret these as slightly overluminous 
tip- AGB stars as there is an accepted physical mechanism 
for doing so: Hot Bottom Burning (HBB). The extra lumi- 
nosity is generated from nuclear burning at the bottom of 
the deep convective zone in massive (^4 M Q ) AGB stars 
(Boothroyd & Sackmann 1992). 

The location of the (M)S-type and carbon stars is as 
expected. Note that the Bertelli et al. (1994) tracks do 
not account for the molecular opacity efects on the atmo- 
spheres of carbon stars which are, however, incorporated 
in more recent models by Marigo (2002). We should also 
emphasize that especially the estimated temperatures of 
the (M)S and C-type stars are uncertain in a systematic 
way, and that in particular the positions of the warmest 
MS-type star (HV 12070) and the warmest carbon star 
(DCMC J050738.74-733233.4) in the diagram may need 
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Fig. 13. The Hcrtzsprung-Russcll diagram for the M-type 
stars (circles), MS or S-type stars (triangles), and carbon 
stars (squares), where the sizes of the symbols are loga- 
rithmically proportional to the mass-loss rate. 



to be corrected towards lower temperatures — otherwise 
the implied Main-Sequence progenitor mass of 5-6 M© 
seems inexplicably large. These two objects have spectra 
that are consistent with a C/O ratio (very) close to unity 
which would cause a reduced opacity in the photosphere. 

Notwithstanding this, it is clear from the diagram (Fig. 
12) that the carbon stars in our sample are descended 
from low-mass Main Sequence stars (mostly in the range 
of roughly 1-2 M©) whereas the M-type AGB stars have 
more massive progenitors (typically 3-5 M Q ). This is as 
expected, as HBB prevents the formation of massive car- 
bon stars (Boothroyd & Sackmann 1993; cf. van Loon et 
al. 2001a; Bergeat, Knapik & Rutily 2002; Marigo, Girardi 
& Chiosi 2003). S-type stars (including stars of MS or CS- 
type) are thought to represent an intermediate phase in 
the formation of a carbon star from an oxygen-rich star 
earlier on the AGB (Lloyd Evans 1983). It is therefore not 
surprising that the (M) S-type stars in our sample are rel- 
atively warm and similarly luminous as the carbon stars 
(i.e. fainter than most M-type stars in our sample). 

The mass-loss rate can be seen to vary across the 
Hcrtzsprung-Russcll diagram due to the evolution of the 
stars and the dependence of the mass-loss rate on lu- 
minosity and stellar effective temperature (Fig. 13). The 
(M) S-type stars and the relatively warm and faint M-type 
AGB stars have the lowest mass-loss rates in this diagram, 
whilst the coolest and most luminous carbon stars and es- 
pecially M-type AGB stars and red supergiants experience 
the most intense mass loss. This illustrates the validity of 
the M(L, T ff) recipe we derived earlier from these data. 
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Fig. 14. Same as Fig. 13, but now the sizes of the sym- 
bols are proportional to (T e ff/Td U st) 2 , which is a measure 
of the distance to the star of the dust formation region, 
-Rdust/ R*- 



There is a clear tendency for the dust shells to be 
warmer for cooler stellar photospheres (Fig. 14). Although 
this may first seem odd, it can probably be understood by 
the dust formation region moving in closer to the star 
as the envelope becomes denser (or vice versa): the dust 
formation radius, i?dust, expressed in stellar radii, i?*, is 
related to the dust temperature at the inner edge of the 
circumstellar dust envelope as Rdust/R* oc (Trff/Tdust) 2 - 
The objects with the relatively coldest dust and largest 
dust-free inner cavity are the early-M type supergiants and 
the coolest and most obscured carbon star in our sample, 
IRAS 05289-6617 (T cff ~ 2200 K and T dust ~ 220 K). The 
latter may have temporarily reduced its mass-loss rate as 
a result of a thermal pulse or (less likely given the short 
timcscale) if it has just become a post- AGB object. 

The fact that the warm and dense circumstellar dust 
envelopes surround the objects with the larger mass-loss 
rates and cooler stellar photospheric temperatures is not 
surprising, as both the mass-loss rate and the cooling of 
the stellar photosphere result from the elevation of the at- 
mospheric layers of the star due to the strong radial pulsa- 
tions which develop on the thermal pulsing AGB and dur- 
ing the red supergiant stage. This may explain the remark- 
able strength of the molecular bands seen in the optical 
spectra of dust-enshrouded M-type AGB stars as well as 
in the optical and near-IR spectra of dust-enshrouded car- 
bon stars in the LMC (van Loon, Zijlstra & Groenewegen 
1999; Matsuura et al. 2002, 2005). 



6. Summary 

We presented new optical spectra for a sample of red 
supergiants and both oxygen-rich and carbon-rich dust- 
enshrouded AGB stars in the LMC. Some of these AGB 
stars have extremely cool stellar photospheres, with spec- 
tral types as late as M10, despite their low metallicity. This 
is cooler than expected on the basis of stellar evolutionary 
models. The same is found for the red supergiants. 

We complemented these new temperature determina- 
tions with similar LMC objects from the literature for 
which spectral subclassifications are available, as well as 
with infrared photometry to sample their spectral energy 
distributions that were then modelled with the radiative 
transfer code dusty to measure the mass-loss rates. The 
mass-loss rates were found to increase with increasing 
luminosity and decreasing stellar effective temperature. 
For the oxygen-rich, M-type AGB stars and red super- 
giants these dependencies were quantified by the mass-loss 
rate formula logM = -5.65 + 1.05 log(L/10, 000 L ) - 
6.31og(T o fl'/3500K), with no evidence for a difference be- 
tween the red supergiants and oxygen-rich AGB stars. 
This formula prescribes the mass-loss rate to within a fac- 
tor two for oxygen-rich stars near the tip of their AGB 
evolution or during the red supergiant phase when they ex- 
perience strong and regular pulsation and enshroud them- 
selves in dusty envelopes, and is valid in particular for the 
superwind phase of extreme mass loss. It thus comple- 
ments the large body of existing formulae derived for less 
massive and less evolved red giants, and can be used in 
models for stellar and galactic evolution in stead of extra- 
polations of those formulae. 

Application of our formula to samples of galactic 
oxygen-rich red giants shows excellent agreement between 
our predictions and mass-loss rates derived from 60 /im 
flux densities for dust-enshrouded AGB stars and the most 
extreme red supergiants. This suggests that the formula 
is valid for stars with masses from Mzams ~1.5M q to 
A^zams >20M Q , and that the total (gas+dust) mass-loss 
rate does not depend sensitively on the initial metal abun- 
dance. On the other hand, our formula predicts higher 
rates than the observational estimates for less extreme red 
giants. This failure may indicate that our formula should 
only be applied to dust- enshrouded objects where the dust 
condensation efficiency is maximal, but it may also be 
explained by errors in the observational estimates if the 
dust-to-gas ratio in these less extreme objects is lower than 
what is commonly assumed. 

Carbon stars and stars of type MS or S seem to show 
a similar behaviour of the mass-loss rate as a function of 
luminosity and temperature, but a quantification of these 
dependencies awaits a more careful and extensive analysis 
of these classes of objects. 

We presented Hertzsprung-Russell diagrams, which 
not only illustrate the trend of mass-loss rate with lumi- 
nosity and temperature but also a trend for dust to form at 
relatively greater distances from early-M type stars than 
in the case of stars with cooler stellar photospheres. These 
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diagrams also showed evidence for Hot Bottom Burning 
preventing carbon star formation for massive AGB stars 
and enhancing their luminosity at the tip of the AGB. 
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Appendix A: Spectra of galactic M-type stars 

During the DFOSC 1995 and EMMI 2000 runs, spectra 
were taken of several galactic M-type stars for reference 
purposes. The observations were performed and reduced 
in exactly the same way as the LMC targets, but with ex- 
posure times ranging between only 5 seconds and 5 min- 
utes. The spectral classification and IR photometry from 
2MASS and IRAS arc listed in Table Al, and the spec- 
tra are displayed in Figs. Al & A2. RSVel was also cov- 
ered by MSX, yielding flux densities of F s .28 = 266.1, 
F 12 .i3 = 244.6, ^14.65 = 188.3 and F 21 . 34 = 121.5 Jy. 
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Table A.l. Spectroscopic targets in the Milky Way, listed in order of increasing Right Ascension (J2000 coordinates) 
and accompanied by the observing run code (cither "D5" = DFOSC 1995, or "EO" = EMMI 2000). Also listed are 
their spectral classifications and near-IR magnitudes and mid-IR flux densities (in Jy, with the wavelengths in /jm). 
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